Cryptographic system for performing secure iterative computations and signal processing directly on encrypted data in untrusted environments

ABSTRACT

Disclosed embodiments include a cryptographic system implemented in at least one digital computer with one or more processors or hardware such as FPGAs for performing iterative secure computations, analysis, and signal processing directly on encrypted data in untrusted environments. According to a basic embodiment, the proposed cryptographic system comprises: (a) at least one secure protocol for performing matrix multiplications in the encrypted domain, and (b) at least one secure iterative protocol for solving systems of linear equations in the encrypted domain. According to a particular embodiment the system comprises a plurality of privacy-preserving protocols for solving systems of linear equations (SLE) directly based on homomorphic computation and secret sharing. More specifically, according to a particular embodiment the system uses a protocol whereby systems of linear equations are solved securely and iteratively without imposing any restrictions on the matrix coefficients.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims the benefit of U.S. Provisional Application No. 61/240,179 filed on Sep. 4, 2009 by the present inventors, which is incorporated herein by reference.

TECHNICAL FIELD

Disclosed embodiments relate to cryptography methods and systems. Specifically, they relate to cryptosystems for the execution of operations directly on encrypted data.

BACKGROUND

In modern society, digital data about individuals can be found relatively easily in the communication networks, especially the Internet. Although the public supports the last decades' advances in digital networks, the sensitive nature of this data motivates the raise of an increasing concern about the public availability of personal data, and the processing performed on it. For instance, in Europe this concern has been reflected in a series of Directives, dealing with the protection of individuals' personal data. Directive 95/46/EC deals with the protection of individuals with regard to the processing of personal data and the free movement of such data, where personal data means any information relating to an identified or identifiable natural person. One of the main ideas of this Directive is that data processing systems must respect the fundamental rights and freedoms, especially the right to privacy. Similar laws and directives regarding informational privacy have been proposed in the US and other developed countries. These include Internet privacy, medical privacy, financial privacy, law enforcement privacy, and political privacy.

There are currently numerous methods available to provide security to a database containing data intended to be held private to third parties. This problem is solved by methods of secure authentication. These methods are widely used to protect the security of databases containing medical, financial, and other types of data. An additional level of security can be added by encrypting all the data contained within the secure database and transmitted to and from the secure database during uploading or downloading operations. This is especially useful when data must be transmitted over a potentially unsecure medium such as computer networks. However, even in these increased security systems at some point the encrypted data is decrypted by an authorized user to perform the necessary data analysis tasks. For instance, in the case of medical research the secure clinical data management system containing the encrypted clinical data and physiological time-series corresponding to a particular study, research, or clinical trial requires the authorized researcher to decrypt the data in order to perform the needed mathematical analysis, signal processing, and statistical analysis in order to generate the study results. There are situations where this conventional framework does not meet the security requirements, that is, in order to analyze and process the encrypted data, such data should not be first decrypted (i.e. the data should be kept encrypted at all times). There are situations where data should be completely private even to researchers, clinical administrators, and system administrators. This may be due to purely privacy reasons or due to research considerations. For instance, certain studies may require scientists and researchers to be completely blinded and the hypothesis and analysis methods to be chosen a-priori. In order to accomplish this, methods for performing the typical mathematical operations and computations including statistical analysis techniques, algebraic methods, signal processing methods, and other computation operations directly on the encrypted datasets are needed. Currently, the availability of such secure methods and cryptosystems is limited.

Conventional cryptographic protocols deal with the problem of protecting some private information from an unauthorized third party that otherwise could modify or have access to the information. In the scenario of secure processing, where the privacy must be preserved not only against an unauthorized third party, but also against the parties that process the data, there are no available systems or methods that can be used in real-world scenarios in terms of both computational cost and communication complexity.

Up to now, the efficient protocols presented in the field of signal processing in the encrypted domain have been focused on linear operations, like scalar products, and non-iterative algorithms. Nevertheless, there are many basic algorithms needed for most signal processing applications that are iterative and involve not only scalar products with known values, but also products between two a-priori unknown sequences. The lack of these algorithms would suppose missing a powerful and irreplaceable tool that enables almost any signal processing application and most types of analysis methods.

Currently there is no practical fully homomorphic cryptosystem, that is, there is no secure cryptosystem that allows for the homomorphic computation of additions and products without restrictions. There has been a recent contribution by Gentry that presents a cryptosystem based on ideal lattices with bootstrappable decryption, and it has been shown that it achieves a full homomorphism. Nevertheless, the authors of this method concede that making the scheme practical remains an open problem.

At the moment there are no cryptosystems available that are capable of performing computations such as solving systems of linear equations without imposing any restrictions on the matrix coefficients. These cryptosystems are a critical building block needed to develop and implement more complex cryptosystems capable of performing advanced signal processing, computation, and analysis directly on encrypted data.

SUMMARY

According to one embodiment, the proposed system involves a cryptographic system (i.e. a cryptosystem) for the execution of operations directly on encrypted data, that is, a system designed to address the problem of efficiently processing signals in untrusted environments, where not only the communication channel between parties is unsecure, but also the parties that perform the computation cannot be trusted.

Disclosed embodiments include a cryptographic system implemented in at least one digital computer with one or more processors or hardware such as FPGAs for performing iterative secure computations, analysis, and signal processing directly on encrypted data in untrusted environments. According to a basic embodiment, the proposed cryptographic system comprises: (a) at least one secure protocol for performing matrix multiplications in the encrypted domain, and (b) at least one secure iterative protocol for solving systems of linear equations in the encrypted domain. According to a particular embodiment, and without limitation, the system comprises a plurality of privacy-preserving protocols for solving systems of linear equations (SLE) directly based on homomorphic computation and secret sharing. More specifically, according to a particular embodiment, and without limitation, the system uses a protocol whereby systems of linear equations are solved securely and iteratively without imposing any restrictions on the matrix coefficients.

BRIEF DESCRIPTION OF THE DRAWINGS

Disclosed embodiments are illustrated by way of example, and not by way of limitation, in the figures of the accompanying drawings.

FIG. 1 shows a block diagram to illustrate the cryptographic system according to one embodiment.

FIG. 2 shows a high level block diagram to illustrate the operation of the cryptographic system.

FIG. 3 shows a block diagram to illustrate the secure iterative protocol for solving systems of linear equations in the encrypted domain according to one embodiment.

DETAILED DESCRIPTION

FIG. 1 shows a block diagram of the cryptographic system 102 according to one embodiment. The embodiments disclosed describe how to realize a cryptographic system 102 for performing secure mathematical computations, signal processing, statistical signal processing, statistics, and data analysis directly on encrypted data 100 in untrusted environments. According to one embodiment, the system involves a cryptosystem 102 for the execution of operations directly on encrypted data 100, that is, a system designed to address the problem of efficiently processing signals in untrusted environments, where not only the communication channel between parties is unsecure, but also the parties that perform the computation cannot be trusted.

Certain specific details are set forth in the following description and figures to provide a thorough understanding of various embodiments disclosed. Certain well-known details often associated with computing and software technology are not set forth in the following disclosure to avoid unnecessarily obscuring the various disclosed embodiments. Further, those of ordinary skill in the relevant art will understand that they can practice other embodiments without one or more of the details described below. Aspects of the disclosed embodiments may be implemented in the general context of computer-executable instructions, such as program modules, being executed by a computer, computer server, or device containing a processor. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. Aspects of the disclosed embodiments may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote storage media including memory storage devices. Those skilled in the art will appreciate that, given the description of the modules comprising the disclosed embodiments provided in this specification, it is a routine matter to provide working systems which will work on a variety of known and commonly available technologies capable of incorporating the features described herein.

According to a disclosed embodiment, the cryptographic system 102 is implemented in a digital computer with one or more processors for performing secure computations and signal processing directly on encrypted data 100 in untrusted environments, said cryptographic system 102 comprises: (a) at least one secure protocol for performing matrix multiplications in the encrypted domain 104, and (b) at least one secure iterative protocol for solving systems of linear equations in the encrypted domain 106. According to a particular embodiment, and without limitation, the system uses a secure iterative protocol whereby systems of linear equations are solved securely and iteratively to directly generate the results 112 of such secure mathematical computations, signal processing, statistical signal processing, statistics, and data analysis directly on encrypted data 100 in untrusted environments. According to other embodiments the disclosed cryptosystem can be implemented in other hardware besides a digital computer including microcontrollers, DSPs, FPGAs or ASICs, as well as in firmware and software.

The disclosed embodiments of the cryptosystem 102 and associated secure protocols have specific and substantial utility by themselves as systems and methods for cryptography and secure (encrypted) computation but also in a variety of practical applications in diverse fields including secure clinical data management systems, secure web-enabled platforms for collaboration involving privacy data, secure law enforcement systems, secure financial systems, and secure military systems involving transmission of encrypted data to be processed in real-time without decryption.

According to a basic embodiment, the proposed cryptographic system 102 comprises: (a) at least one secure protocol for performing matrix multiplications in the encrypted domain 104, and (b) at least one secure iterative protocol for solving systems of linear equations in the encrypted domain 106 whereby systems of linear equations are solved securely using a secure multiplication protocol based on homomorphic computation 108 and secret sharing 110 without placing any restrictions on the nature of the matrix coefficients (i.e. the matrix coefficients can be real numbers or complex numbers).

By way of example, and not by way of limitation, the embodiments of the disclosed cryptosystem 102 can be used to perform direct analysis and processing in encrypted clinical, financial, national security, military, law enforcement, and political data without the need to decrypt said data prior to performing said analysis and processing. For instance, according to one embodiment the disclosed cryptosystem 102 can be used to realize and implement a secure clinical system wherein all patient data, clinical history, clinical data, physiologic data and time series, biochemical data, and drug therapy data are encrypted at all times for secure transmission, management, and collaboration among researchers and clinicians; and all operations, processing, and statistical analysis performed on such data is conducted on the encrypted domain (i.e. with the data encrypted at all times). A specific particular embodiment of such secure clinical data management and analysis system comprises (a) at least one secure protocol for performing matrix multiplications in the encrypted domain 104, and (b) at least one secure iterative protocol for solving systems of linear equations in the encrypted domain 106 in addition to the well-known elements comprising clinical management systems.

FIG. 2 shows a high level diagram to illustrate the operation of the cryptographic system according to one embodiment. In the embodiment illustrated, and without limitation, there are two parties

200 and

220. The first party owns an encrypted matrix [[A]] and an encrypted vector [[b]], as well as the key for producing encryptions using a homomorphic encryption 202. The party

owns both the encryption and decryption key 216 for the same homomorphic encryption. The cryptographic system is used at both ends 206 and 212 in order to privately solve the linear system of equations A·x=b. This cryptographic system comprises a protocol of communication between both ends 208 and 214 that defines the operations performed and the format of the interchanged numbers in order to interactively solve the system without disclosing any information. This communication takes place over any communication channel 210 (e.g., a wired or wireless medium or a sole device). After the execution of the protocol, the system outputs the solution x 204 and 218 to the linear system of equations A·x=b. It is important to note that both

and

can represent a plurality of parties or entities (not necessarily a single entity).

Secure matrix multiplication (linear transformations) and solving linear systems of equations iteratively in the encrypted domain are fundamental mathematical building blocks to solve a great variety of problems involving data analysis, time-series analysis, digital signal processing, statistical signal processing, optimal filtering, adaptive filtering, digital communications, coding, encryption, information theory, and any other problem involving the least squares framework and representation of signals in vector spaces.

While particular embodiments are described, it is understood that, after learning the teachings contained in this disclosure, modifications and generalizations will be apparent to those skilled in the art without departing from the spirit of the disclosed embodiments.

A. Mathematical Notation.

In the detailed description the following mathematical notation is used. The mathematical notation is thoroughly described in this section to aid those skilled in the art to understand, implement, and practice the disclosed embodiments.

In the description we use indistinctly lowercase letters to represent classes in a ring (

+, ·) and a representative of that class in the interval [0, n). ┌·┘ represents the rounding function of a number to the nearest integer.

The used vectors have size L and are represented by lower-case boldface letters, whereas matrices are represented by upper-case boldface letters. A′={a_(i,j)}_(r,s) ^(t,u) represents the submatrix of A of size (t−r+1)×(u−s+1), defined by a′_(i,j)=a_(i+r−1,j+s−1).

The encryption of a number x is represented by [[x]], and the vector (matrix) formed by the encryptions of the vector x (matrix X) is represented by [[x]] ([[X]]).

The operations performed between encrypted and clear numbers are indicated as if they were performed in the clear; e.g. [[X]]·b represents the encryption of [[X·b]].

Regarding the complexity calculations, the complexity of basic modular operations, like additions (A), products (P) and exponentiations (X) is denoted by Comp_(A), Comp_(P), Comp_(X) respectively, prefixing an E (i.e. EA, EP, EX) when they are performed under encryption. The factor ct<1 denotes the ratio between the size of a clear-text value and that of an encrypted value. Finally, the subscript cm denotes communication complexity, measured in number of sent encryptions, while cp indicates computational complexity, with an indication of the party whose complexity is represented.

B. Techniques Employed According to One Embodiment.

According to one embodiment, the cryptographic system 102 uses secure multiparty computation techniques including homomorphic encryption and secret sharing.

According to one embodiment the cryptographic system 102 makes use of homomorphisms between the groups of clear-text and cipher-text, that allow for the execution of a given operation directly on encrypted values, without the need for decryption. This includes homomorphic RSA cryptosystem, with a multiplicative homomorphism, or Paillier with an additive homomorphism. The methods implemented are not restricted to the use of cryptosystem for the presented protocols, as far as it presents an additive homomorphism. In this embodiment the system 102 uses an extension of Paillier encryption in its threshold and non-threshold form; that is, a k out of M threshold public key encryption system is a cryptosystem where the private key is distributed among M parties, and at least k of them are needed for decryption.

According to one embodiment the cryptographic system 102 makes use of secret sharing. In this embodiment a given value (the secret) is divided among several parties, such that the cooperation among a number of these parties is needed in order to recover the secret. None of the parties alone can have access to the secret. According to one embodiment the scheme is based on polynomials, and the need of k points in order to completely determine a degree (k−1) polynomial. The disclosed embodiment uses secret sharing for two-party protocols based on linear functions, without limitation, since each party can represent a plurality of entities. According to one embodiment, the system supports the computation of sums and products directly on the shares as follows: 1) let

be the domain of the secrets; 2) then, a share of a secret x is defined as two values x_(A) and x_(B), owned by their respective parties, such that x_(A)+x_(B)≡x mod n; 3) hereinafter, randomizing an encrypted value x means obtaining one share and providing the encryption of the other (through homomorphic addition).

C. Definitions of Protocols and Constructions According to One Embodiment.

According to one embodiment of the cryptosystem 102, and without limitation, the protocols are based on two parties,

and

both using an additively homomorphic cryptosystem in an asymmetric scenario, where

can only encrypt, but

also possesses the decryption key, and can perform both encryption and decryption. For the problem of solving an SLE A·x=b, the system requires that

owns an encrypted version of the system matrix [[A]], and of the independent vector [[b]].

According to one embodiment, and without limitation, the SLE may have A being either a positive definite matrix or a strictly diagonally dominant matrix. In this particular case, it can be guaranteed both a solution to the system and the convergence of the studied methods, as detailed later.

Regarding the privacy requirements, according to one embodiment the cryptosystem 102 is especially adapted for semi-honest parties, in the sense that they adhere to the established protocol, but they can be curious about the information they can get from the interaction. In this embodiment, the protocols can be proven private; informally, both parties

and

can only get the information leaked from the solution to the system, and no information is leaked from the intermediate steps of the protocols.

D. Secure Computational Engine and Method for Multiplying Matrices in the Encrypted Domain.

In order to multiply two encrypted matrices, as there is no multiplication operation in an additively homomorphic cryptosystem, it is necessary to execute an interactive protocol in order to perform each product. According to one embodiment, the systems uses the following method for non-threshold encryption.

According to one embodiment of the secure protocol, and without limitation, it is assumed that there exists an additively homomorphic cryptosystem with plaintext in

such that

can decrypt and both

and

can encrypt. In this embodiment,

owns two encrypted scalars [[x₁]] and [[x₂]] and wants to multiply them under encryption. In order to do that, according to one embodiment

generates two random values r₁, r₂ε

and uses them to blind both numbers, using homomorphic modulo-n sum obtaining [[z₁]]=[[x₁]]+r₁ mod n, and [[z]]=[[x₂]]+r₂ mod n, and sends them to

.

According to one embodiment of the method, due to his decryption capabilities,

can obtain z₁ and z₂ in the clear, multiply them, and reencrypt the result [[z₁·z₂]].

sends this encrypted product to

who, through homomorphic sums, can obtain the desired result, as [[x ₁ ·x ₂ ]]=[[z ₁ ·z ₂ ]]−r ₁ [[x ₂ ]]−r ₂ [[x ₁ ]]−r ₁ r ₂.

According to the embodiment involving a threshold homomorphic cryptosystem, the procedure is analogous, with the exception that the random values must be generated by both parties.

According to the embodiment involving a product of an L×M matrix and a scalar, the protocol is exactly the same as the scalar-scalar case, with L×M scalar products in parallel.

According to the embodiment involving a matrix-matrix product, all the scalar products are performed using the scalar-scalar product protocol in parallel, with only one randomization per matrix coefficient, and the remaining operations are sums, that can be performed homomorphically. According to one embodiment, in order to minimize the computation and communication complexity,

may let

perform all the partial additions that

can do in the clear and

would need to do homomorphically.

Neglecting the complexity of the random number generation algorithms, the complexity of the whole protocol, when multiplying an L×M matrix and an M×N matrix is Comp_(cmMULT)(L,M,N)=M·(L+N)+L·N Comp_(cpMULT,A)(L,M,N)=L·N·M·(3Comp_(EA)+2Comp_(EP)) Comp_(cpMULT,B)(L,M,N)=M·(L+N)Comp_(Decrypt) +M·L·N·NComp_(P) +L·N·((M−1)Comp_(A)+Comp_(Encrypt)). E. Secure Computational Engine and Iterative Protocol for Solving Linear Equations in the Encrypted Domain.

FIG. 3 shows a block diagram of an illustrative method for solving systems of linear equations iteratively in the encrypted domain according to one embodiment. According to one embodiment, the cryptosystem 102 includes an iterative method (protocol) for solving systems of linear equations directly on the encrypted data 100 (i.e. without the need to decrypt the data). According to one embodiment presented here for illustrative purposes only, and not by way of limitation, the cryptosystem includes an iterative method for solving systems of linear equations as described as follows. Let L denote the dimension of the system matrix A.

blinds the encryption of the principal diagonal of [[A]] through homomorphic addition 300.

sends these encryptions to

302, so that

can decrypt a share of the principal diagonal of [[A]] 304. Then,

and

calculate, through parallel secure multiplication protocols, the encrypted diagonal matrix [[γD⁻¹]] and the factor γ 306.

and

using the secure multiplication protocol, can calculate encryptions of the modified matrix [[γM]] and the modified vector [[γc]] 308. After that,

and

end up with shares of [[γM]] 310.

performs the first iteration of the protocol using homomorphic operations 312, obtaining [[x′₁]]=[[δx⁽¹⁾]]. Then, for each iteration, being k the number of the current iteration, both parties follow the next procedure:

and

execute the secure product protocol on [[δM]] and [[x′_(k)]] 314; using the result of the previous operation,

homomorphically calculates the encryption of [[x′_(k+1)]] 316. After each iteration or after a predefined number of iterations, a stopping condition is checked 318. When the stopping condition 318 is fulfilled, the protocol stops, and

sends the encrypted [[x′_(k)]] to

for decryption 320, and both parties obtain the resulting x^((k)) 322.

According to one embodiment, the general form of stationary iterative methods for solving SLEs is x ^((k+1)) =M·x ^((k)) +c.

According to one embodiment, and without limitation, the system matrix is decomposed into A=D(L+I+U), a diagonal matrix D, a lower triangular matrix L and an upper triangular matrix U, having both L and U zeros in their principal diagonals. Then, M=−(L+U) and c=D⁻¹b. We refer to this embodiment as the Jacobi method, Jacobi embodiment, Jacobi protocol, or Jacobi secure protocol due to the nature of the matrix decomposition and the iterative method for solving the system of linear equations. As divisions are not supported homomorphically, the previous iteration cannot be implemented directly. Thus, according to this embodiment, the division is simulated by multiplying each row of A by the diagonal elements of the remaining rows, what results in multiplying the matrix M of Jacobi method by a scalar factor

$\gamma = {\left( {\prod\limits_{i = 1}^{L}\; a_{ii}} \right).}$ A′=γD ⁻¹·(A−D)=γM.

According to the Jacobi embodiment, the factor γ is propagated at every iteration of the algorithm: γ^(k) x ^((k)) =−γD ⁻¹(A−D)·γ^(k−1) x ^((k−1))+γ^(k) D ⁻¹ b.

According to one embodiment the secure iterative protocol assumes that

can decrypt and both

and

can encrypt with an additive homomorphic scheme, and that

owns encryptions of [[A]] and [[b]] with this homomorphic system. In order to allow for efficient computation, the following protocol is executed:

-   -   1.         can blind the principal diagonal of [[A]] and send it to         .     -   2.         decrypts it, both parties ending up with additive shares of the         diagonal elements {a_(ii)}.     -   3. With this shares, both parties can securely compute shares of         the diagonal matrix

${\left( {\gamma\; D^{- 1}} \right)_{jj} = {\prod\limits_{\underset{i \neq j}{i = 1}}^{L}a_{ii}}},$

-   -    through ┌ log₂(L−1)┐ rounds of parallel secure multiplication         protocols. They can also calculate the value of γ, and disclose         it for use in the following steps of the protocol.     -   4.         can then calculate the encryption of [[γM]]=[[−γD⁻¹]]·[[A−D]]         and [[γc]]=[[γD⁻¹]]·[[b]], invoking the secure multiplication         protocol.     -   5. Then,         sends         a blinded and encrypted version of [[γM]], that         decrypts for use in the following iterations.

After these initial steps, for the first iteration of the secure protocol both parties agree in an initial vector x⁽⁰⁾ and

calculates, through homomorphic additions and multiplications, the encryption of [[γx⁽¹⁾]]=[[γM]]·x⁽⁰⁾+[[γc]].

For each subsequent iteration,

calculates the encryption of γ·[[γ^(k−1)c]], and then both parties use the secure multiplication protocol and homomorphic additions in order to obtain the vector for the following step [[x′ _(k)]]=[[γ^(k) x ^((k)) ]]=[[γM]]·[[γ ^(k−1) x ^((k−1))]]+γ·[[γ^(k−1) c]].

It must be noted that the matrix [[γM]] does not have to be communicated at each iteration, as its blinded version was stored by B at the initial step. Thus, only two vectors per iteration are sent between

and

.

After each iteration, the factor γ multiplies the result; thus, after a number of steps, the cipher will not be able to accommodate the scaled number, and the protocol will have to stop. It must be noted that the accumulated factor is not only γ, but also the quantization step Δ used for the initial quantization of the coefficients of both the system matrix A and the vector b in order to make them integers so that they can be encrypted. This factor must also be taken into account every time γ multiplies vector c, so that the homomorphically added vectors be quantized with the same scaling factor.

According to one embodiment, each step of the protocol according to this embodiment can be proven secure with semihonest parties, due to the semantic security of the underlying cryptosystem, the security of the multiplication protocol, and the fact that the unencrypted values that each party sees are random and uncorrelated.

The complexity of the initial part (Jacobi Initial, JI) of the protocol is Comp_(cmJI)=3L ²+2┌ log₂(L−1)┐−3L+5+ct Comp_(cpJI,A)=(5L ²+4┌ log₂(L−1)┐−5L+8)Comp_(EA)+(L ² +L┌ log₂(L−1)┐−L+2)2Comp_(EP) Comp_(cpJI,B)=(L ² +L┌ log₂(L−1)┐−L+2)(Comp_(Decrypt)+Comp_(P)+Comp_(Encrypt))+(L ² −L+1)Comp_(Decrypt).

The first iteration (J1) does not involve any interaction, and

incurs in a computational complexity of Comp_(cpJ1,A)=L²(Comp_(EP)+Comp_(EA)).

The complexity of each of the subsequent iterations of this protocol (J) is the following Comp_(cmJ)=2L Comp_(cpJ,A)=(3L ²−2L)Comp_(EA)+(2L ² −L)Comp_(EP) +LComp_(EA) Comp_(cpJ,B) =L(Comp_(Decrypt)+Comp_(Encrypt))+(L ² −L)Comp_(P)+(L ²−2L)Comp_(A).

After a number of iterations, either the solution can be disclosed, or an error metric can be obtained to determine whether convergence has been achieved. While the choice of this error metric is arbitrary, one possibility according to one embodiment is to homomorphically subtract [[x^((k))]]−[[x^((k−1))]], and either decrypt the result or perform L parallel encrypted comparisons with a predetermined threshold.

According to one embodiment, and without limitation, the coefficients of the system matrix A are quantized versions of the real-valued coefficients, with a quantization step Δ, such that their quantized absolute value is bounded by an integer T>0.

For the first part of the iterative protocol embodiment disclosed, where the factor γ and the matrix γD⁻¹ are calculated, γ is the highest number that the system will have to represent, and it is bounded by T^(L); the bound for the elements of γD⁻¹ is T^(L−1). Furthermore, as γ is disclosed in the following step, it can constitute a more accurate bound to the encrypted coefficients of γD⁻¹. A bound for the absolute value of the coefficients of γM and of γc is min(T^(L), γT).

According to one particular embodiment, the bound to which the elements of the first vector x⁽¹⁾ are subject is L·T^(L+1).

In each iteration, the previous bound is multiplied by L·T^(L), meaning that the bound for the elements of the k-th iteration is L^(k)·T^(k·L+1), i.e., the needed bit-size of the cipher is linear both in the dimension of the system and in the maximum number of iterations that can be performed without errors. Furthermore, the quantization step of the elements of x^((k)) will be Δ^(kL).

When dealing with iterative algorithms such as the embodiment presented here, it is necessary to determine whether the algorithm can converge or not before applying the algorithm. In the general case of stationary iterative methods, the necessary and sufficient condition for their convergence with an arbitrary initial vector x⁽⁰⁾ is that max_(i)|γ_(i)(M)|<1, where λ_(i)(M) are the eigenvalues of M. In this particular embodiment, M=−D⁻¹·(A−D). Let us assume that A is a strictly diagonally dominant matrix with bounded coefficients |a_(ij)|≦T. By Ostrowski's theorem, the eigenvalues of M are located in the union of L discs

${\mathcal{L}_{1}\overset{\Delta}{=}{\bigcup\limits_{i = 1}^{L}\left\{ {z \in {{\mathbb{C}}:{{{z - m_{ii}}} \leq {\min\left\{ {R_{i},C_{i}} \right\}}}}} \right\}}},{where}$ ${m_{ii} = 0},{m_{ij} = \frac{a_{ij}}{a_{ii}}},{i \neq j},{and}$ ${R_{i} = {\sum\limits_{{j = 1},{j \neq i}}^{L}{m_{ij}}}},{C_{i} = {\sum\limits_{{j = 1},{j \neq i}}^{L}{{m_{ji}}.}}}$

As A is strictly diagonally dominant, Σ_(j=1,j≠i) ^(L)|a_(ij)|<|a_(ii) 651

R_(i)<1. Thus, it is possible to bound the moduli of the eigenvalues of M as |λ_(i)(M)|<1.

Consequently, the Jacobi embodiment disclosed always converges for strictly diagonally dominant matrices, and the test of convergence is not needed.

F. Extension to Systems with Complex Elements According to One Embodiment.

The protocol presented in the earlier section does not impose any restrictions on the matrix coefficients. It can be used for both unrestricted real coefficients and complex coefficients. In this section we explain how the protocol can be implemented in the case of complex coefficients.

For convenience, and without loss of generality, a binomial representation is chosen for complex numbers; the encryption of a complex number [[x]] is given by the pair [[x]]

([[

{x}]], [[

{x}]]), that is, the pair of encryptions of its real and imaginary parts.

According to one embodiment, the complex addition operation can be performed through two real additions; thus, it can be performed homomorphically between two encrypted complex numbers as two homomorphic real additions. Regarding the complex product, when it involves one known factor and one encrypted number, it can also be performed homomorphically as four homomorphic multiplications and two homomorphic additions. When involving two encrypted factors, it is performed through four parallel (real) scalar multiplication protocols and two homomorphic additions.

Thus, the matrix multiplication protocol disclosed in the previous section works on complex numbers by adopting the proposed encrypted complex representation, and by substituting the real additions and products by their corresponding complex operations. As a result, the communication complexity is doubled, the number of performed products is multiplied by four, and the number of performed additions is doubled.

For the secure SLE solving protocol, the needed modifications are essentially the same, being the involved scale factors also complex. We must remark that the hypotheses on the system matrix that were imposed for ensuring the convergence of the algorithms are kept unaltered when dealing with systems with complex coefficients, as they are general hypotheses not restricted to systems with real coefficients.

While particular embodiments have been described, it is understood that, after learning the teachings contained in this disclosure, modifications and generalizations will be apparent to those skilled in the art without departing from the spirit of the disclosed embodiments. It is noted that the foregoing examples have been provided merely for the purpose of explanation and are in no way to be construed as limiting. While the system has been described with reference to various embodiments, it is understood that the words which have been used herein are words of description and illustration, rather than words of limitations. Further, although the system has been described herein with reference to particular means, materials and embodiments, the actual embodiments are not intended to be limited to the particulars disclosed herein; rather, the system extends to all functionally equivalent structures, methods and uses, such as are within the scope of the appended claims. Those skilled in the art, having the benefit of the teachings of this specification, may effect numerous modifications thereto and changes may be made without departing from the scope and spirit of the disclosed embodiments in its aspects. 

The invention claimed is:
 1. A homomorphic cryptographic system, comprising: (a) a memory configured to store computer-executable instructions; (b) at least one hardware processor configured to execute said computer-executable instructions that when executed cause the processor to: (1) perform matrix multiplications in an encrypted domain using a multiplication protocol based on homomorphic computation and secret sharing, wherein data remains encrypted at all times within the encrypted domain by: (i) generating a set of random values; and (ii) obtaining an encrypted product of encrypted matrices using said set of random values and homomorphic operations; (2) solve systems of linear equations in the encrypted domain using at least one secure iterative protocol configured for operating with matrices containing no restrictions on the coefficients being real or complex by: (i) generating shares of a diagonal matrix through homomorphic operations; (ii) calculating an encrypted modified matrix and an encrypted modified vector using secure multiplication; and (iii) performing iterations using homomorphic operations; wherein said homomorphic cryptographic system performs iterative computations, signal processing, and data analysis on encrypted data in untrusted environments without the need to decrypt said encrypted data.
 2. The homomorphic cryptographic system of claim 1, wherein said secure iterative protocol for solving systems of linear equations in the encrypted domain is based on privacy-preserving protocols for homomorphic computation and said iterative protocol is configured for operating with elements of an infinite field including matrices with real or complex coefficients.
 3. The homomorphic cryptographic system of claim 2, wherein said secure iterative protocol for solving systems of linear equations in the encrypted domain is based on privacy-preserving protocols based on homomorphic computation and secret sharing.
 4. The homomorphic cryptographic system of claim 3, wherein said secure iterative protocol for solving systems of linear equations in the encrypted domain is based on privacy-preserving protocols based on homomorphic computation and secret sharing, and said secure iterative protocol for solving systems of linear equations is implemented as a stationary iterative method.
 5. The homomorphic cryptographic system of claim 4, wherein said secure iterative protocol for solving systems of linear equations in the encrypted domain is based on a secure protocol employing an iterative Jacobi update matrix.
 6. The cryptographic system of claim 5, wherein said hardware processor is further configured to perform the steps of: (a) calculating one or more shares of a diagonal of a system matrix; (b) calculating an encrypted modified matrix and an encrypted modified vector; (c) calculating an encrypted output vector based on homomorphic computation and secret sharing; and (d) decrypting one or more elements of said output vector. 